CLNS 04/1898 



|hep-ph/041 10271 

November 1, 2004 



Impact of Four-Quark Shape Functions on 
Inclusive B Decay Spectra 



Matthias Neubert 

Institute for High-Energy Phenomenology 
Newman Laboratory for Elementary-Particle Physics, Cornell University 

Ithaca, NY 14853, U.S.A. 

and 

School of Natural Sciences, Institute for Advanced Study 
Princeton, NJ 08540, U.S.A. 



Abstract 

It has recently been pointed out that a new class of subleading shape functions in- 
volving i?-meson matrix elements of non-local four-quark operators contributes at order 
^QCD/iTT'b to B ^ Xu I'v decay distributions in the endpoint region. The corresponding 
functions fu{^) and fv{^) are estimated using the vacuum-insertion approximation. A 
numerical analysis of various B —>■ Xu l~v decay spectra suggests that these power cor- 
rections are very small, below present theoretical uncertainties due to other subleading 
shape-function contributions. 



Inclusive decays of B mesons into final states containing light particles, such as B — > 1"^ 
and B ^s7, play an important role in the extraction of the element \Vub\ of the quark mixing 
matrix. Experimental cuts in the analysis of these processes restrict the hadronic final state to 
have large energy, Ex ^ rriB, but only moderate invariant mass, Mx ~ ^/^b-^qcb- In this re- 
gion of phase space, the inclusive rates can be calculated using a twist expansion, which resums 
infinite sets of local-operator matrix elements into non-perturbative shape functions mi21- It 
is well known that the leading term in this expansion obeys a QCD factorization formula 
which separates contributions associated with the hard scale rrib, the jet scale y^mbAqcD, and 
the soft scale Aqcd- In recent work, the calculation of next-to- leading perturbative corrections 
to the various components in this formula has been completed |3] . 

Remarkably, for the case of inclusive decay distributions QCD factorization can be extended 
beyond the leading order in the heavy-quark expansion. Using the formalism of soft-coUinear 
effective theory (SCET) [3111, it has been argued that a factorization theorem holds at every 
order in Aqcd/''^6 |H1 HI- Inclusive decay spectra therefore provide an example of a class 
of observables which have a systematic expansion in non-local string operators, built out of 
quark and gluon fields with light-like separation connected by Wilson lines. Non-perturbative 
hadronic physics is encoded in forward 5-meson matrix elements of these operators in heavy- 
quark effective theory (HQET) fO] • The expansion in string operators is a generalization of the 
conventional (local) operator product expansion for correlation functions at large (Euclidean) 
momentum transfer. 

Recently, two groups have presented the first complete analyses of subleading power correc- 
tions to arbitrary B I'v decay distributions in the shape-function region |H|in|- In these 
studies, the QCD current- current correlator, whose forward 5-meson matrix element is related 
to the hadronic tensor W^^ , is matched onto correlation functions in SCET. These correlators 
are then expanded in terms of non-local string operators in HQET. The resulting expressions 
generalize (and in some cases correct) previous results in the literature IT^ IT^ . In par- 
ticular, it was pointed out that there exist some tree-level contributions involving subleading 
shape functions defined in terms of non-local four-quark operators, which had not been con- 
sidered previously. Their contributions to arbitrary decay distributions can be parameterized 
in terms of two functions fui!^) and fv{uj) defined as P 



/ due-'^' \Uu)T^ + f,{uj)T^ = {-if fdt^ f 



^^^ {B{v)\0,,{t,M,t)\B{v)) 
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where 

Oigih, t2,t) = {hS)o Ti i 7^ tA iS^q)t,n {qS)t2n 7± i Tj {S^h)tn (2) 

is a non-local string operator with fields ordered along the light-cone defined by a vector n. 
Here h is a. heavy-quark field in HQET, g is a massless quark {q = u for semileptonic decay), 
S denotes a soft Wilson line, and tA are the generators of color SU(A'"c). The light-like vector 
n points in the direction of the final-state hadronic jet X^- We work in the 5-meson rest 
frame, where f ^ = (1, 0, 0, 0), and take n^^ = (1, 0, 0, 1). Perpendicular Lorentz indices refer to 
the transverse plane orthogonal to v and n (see jS] for further definitions and notation). The 
traces 
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depend on the Dirac structures of the flavor- changing weak currents j/ = 6 Fj g and Jj = qTjb 
in the deflnition of the hadronic tensor, and Fj, F^ are in principle arbitrary Dirac matrices. 
Expanding both sides of (P) in powers of t, one flnds that the functions fu and have zero 
norm and flrst moment, whereas their second moments are given in terms of a local four-quark 
matrix element, 



(mm i if tA q qii i F, tA h\B{v)) 
f^[u})Ti + f\,[u)Ti = ^ — . (4) 



dujuj 

The subleading shape functions /„ and fy enter the decay rates multiplied by a factor nasifJ^i) ~ 
1, where fii ~ A/'^bAqcD is an intermediate matching scale. Because of the factor tt, this 
perturbative coupling does not provide a numerical suppression. 

While the authors of |HJ agree on the structural form of the four-quark contributions, 
they differ in the assessment of the expected numerical importance of their effects. In [U] , it is 
argued that the four-quark contributions are expected to be suppressed with regard to other 
subleading shape functions. On the other hand, the authors of [8J employ arguments based 
on naive dimensional analysis (including counting factors of 4) to speculate that these effects 
may be the dominant source of Aqcd/'"^?) corrections, which could lead to 0(1) effects in some 
decay spectra. In light of this controversy, it may be of some value to have well-motivated, if 
model-dependent estimates of the subleading shape functions fu and f^. This is what we will 
provide in this Letter. 

We flrst consider the case where the flavor of the light quark q matches that of the B- 
meson spectator quark. It is empirically well established that the magnitude of forward 
5-meson matrix elements of local four-quark operators can be estimated by inserting the 
vacuum intermediate state, thereby factorizing them into products of current matrix ele- 
ments. This approximation is routinely used, e.g., in the analysis of lifetime ratios of beauty 
hadrons IT^ . Depending on the color structure of the operators, it is conventional to 
deflne {Bg\hT^qqT2h\B,) = Bi {Bg\hT^ q\0) {0\qT2 h\Bg) and Fi gg Fs = 

Ei {Bg\hTi q\0) (0|gF2 h\Bq), where the subscript "i" refers to different Dirac structures. The 
\a.Tge-Nc counting rules of QCD imply Bi = 0(1) and = 0{1/Nc), and this hierarchy is 
preserved under renormalization. Theoretical work based on lattice QCD ^Hl and QCD 
sum rules fTH^ suggests that the Ei parameters are indeed rather small, typically of order 0.1 
or less. (More speciflcally, [T7] quotes £i ~ £2 ~ 0.01 at a scale of 2.7 GeV, whereas |18J flnds 
El ^ —0.04 and E2 ~ 0.06 at a scale of 1 GeV.) This is consistent with empirical flndings. For 
instance, at a renormalization point fi ~ "^6/2, the lifetime ratio of charged and neutral B 
mesons can be written as ^H] 

t(B+) I 

-j—^ ^ 1 + 0.0445i + 0.00352 - 0.74ei + 0.20^2 = 1-086 ± 0.017 , (5) 

t(-D ) 

where the last result is the current experimental value. If the Ei parameters were much larger 
than of order 0.1, a flne tuning would be required in order to avoid a large deviation from the 
experimental value. On the other hand, with Bi ^ 1 and Ei = 0(0. 1) it is easy to reproduce 
the experimental result. 
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There appears to be no reason why color suppression should be less effective for non-local 
four-quark operators than for local ones. As a model, we are thus led to replace 



{B^{v)\0^,{hM,t)\B,{v)) ^ e {B^{v)\h{Q) [0,ti] r,r/7^^ g(tin)|0) 

X (0|g(t2n) [t2, t] 7^ i T,- h{tn)\B^{v)) , 



(6) 



where [tk^U] denotes a straight soft Wilson line connecting the points tkU and t^n, and e is 
the color-suppression factor. This ansatz completely specifies our model for the subleading 
shape functions. Whereas the definitions of the parameters Bi and Si introduced in [1^] are 
completely general, the introduction of the parameter e in the equation above corresponds to 
a model hypothesis, because we assume that e is independent of the position arguments ti, 
t2, t and of the Dirac structures F,, Tj. Nevertheless, as we shall see, our model provides an 
expression for the subleading shape functions and with all the right properties such as 
the correct support and moment relations, as far as they are known. We expect that it predicts 
the rough overall scale of the effect reliably. In fact, many analyses of inclusive B decays are 
based on measurements of partial decay rates over kinematical domains that fall in between 
the shape-function region and the region where a conventional operator product expansion 
can be applied In such a situation, the dominant contributions to the decay rates are 
associated with the lowest non-zero moments of the shape functions. In the case of the four- 
quark shape functions, the lowest non-zero moments are given by the matrix element of the 
local four-quark operator in for which the vacuum insertion approximation is certainly 
reasonable. 

The B vacuum matrix elements of the non-local quark bilinears in © can be expressed 
in terms of the leading light-cone distribution amplitude of the B meson 19j, 



(0|g(tn) [t' ,t\iT h{tn)\Bq{v)) 



tr 
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where we have used that the B meson in HQET carries momentum Ai;, with A = — m^. 
Throughout this note we work at lowest order in perturbation theory, so that we can ignore 
the scale dependence of the various objects in the above equation. Taking into account that 
the distribution amplitude 0+ is real, we obtain in our model 



/ dt2{B^{v)\Oi^{h,t2,t)\Bg{v)) 



■iAt 



tr 



■ tr 



7±' 
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The product of traces can be expressed in terms of the objects T\ and T4 defined earher, using 



the fact that between 



'1 + y) any Dirac matrix can be decomposed into the basis 1, 75, 



j^^. A straightforward calculation shows that 



Ti + r4 = - tr 



75^4-^riyi7 



tr 



7^ 



■75 



(9) 
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Taking the Fourier transform of (P), we then obtain 



where 



g{uj) = 6iA-uj) 



+ 



A - u; 



K-uo -uj' 



(10) 



(11) 



'P" denotes the principal-value prescription. The function g has support on the half-interval 
-oo < cij < A, as is indeed required for all shape functions P3I2|- Its first few moments are 



doj giuj) = 



duouo g{uo) = 



du uj'^giuj) = 1 . 



(12) 



The first two conditions ensure that the subleading shape functions fu and have vanishing 
norm and first moment, as required on general grounds [H]. 

In order to have some exphcit models at hand, we adopt two forms of the 5-meson light- 
cone distribution amplitude motivated by QCD sum rules fHl 12015 namely 

2 



iGN, 



00 
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UJo 
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— A, 
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(13) 



where the relation between c^o and A is implied by the equations of motion. Introducing the 
dimensionless variable z = (A — uj)/ujq = lo/ujo > 0, we obtain 



1 



[6{z) -2e'' + 2ze-^'Ei{z)] , 

(.ln.+ |(l + .2)-^(l-.2)) 



4 ^ 



[l + z 



(14) 



where Ei(z) = — P j'^^dt /t is the exponential-integral function. The variable Cj = K — uj 
is the most convenient one when calculating decay spectra, because it is independent of the 
definition used for the 6-quark mass [5 . A graphical representation of these two functions is 
shown in Figure H Notice the rapid fall-off of the model functions for values uj away from the 
endpoint. For the first model, this reflects the assumed exponential fall-off of the distribution 
amplitude 0+. But even in the second model, for which has only a power-like fall-off, the 
function g decreases quickly in magnitude. 

So far, we have assumed that the light quark q is contracted with the spectator quark in 
the B meson. But this is not always possible (e.g., if they have different flavor), and one may 
entertain the possibility that the light-quark pair in the operator O^^q is produced from the 
vacuum. In this possible factorization model is to write 

{B{v) I {hS), Fi tA {S^q)t,n {qS)t,n {S^h\r^\B{v)) 

^B{B{v)\[{hS),T,tAY^ [T,tA{S^h),^Y\B{v)) (0| [(5tg),,„];^ [(g5)*.n]>), (15) 
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Figure 1: Model predictions for the function g{ijj) obtained using A = 0.63 GeV, 
corresponding to = 4.65 GeV. The thick hues represent the 5-function locahzed 
at = A. 



where i3 is a bag parameter, and we have left the Dirac structures arbitrary for the moment. 
Lorentz and gauge invariance dictate that 



(0|[(^^gkn]:[(g^)*.n];|0) 



{qq) 1 + iAhit2 - ti)i 



(16) 



where (qq) is the local quark condensate, and is a hadronic parameter. From above, we 
now obtain in our model 

{B{v) I ChS)o Fi tA {Sk)t,n {qS)t,n ^2 tA {S^hU\B{v)) 



{qq) {B{v)\h{Q)[Q,t]V,T2h{tn)\B{v)) 



AN. 



+ iKh{t2 - h) {B{v)\h{0) [0,t] Fi ^F2 h{tn)\B{v)) 



;i7) 



Applying this formula to the specific case of the operator in (Q), for which Fi = Fj]/i7^ and 
F2 = 7^fiFj, we see that the Dirac structures Fi F2 and Fi yiF2 both vanish due to yi^ = 0. 
Therefore, the non-valence contributions vanish in the vacuum-insertion approximation. 

Let us briefly discuss the phenomenological implications of our results by considering three 
decay spectra in semileptonic B —>■ V„ l~u decay, referring to jH] for details and derivations. In 
all cases we show the contributions of the leading-order shape function S{i2}) and of the four- 
quark shape functions as calculated in our model, ignoring other subleading shape-function 
contributions, which have already been discussed in [Hj. In our model, only decay distributions 
of charged B mesons are effected by the four-quark contributions. Of particular interest for a 
measurement of \Vub\ are the spectra in the variables Ei (charged- lepton energy), P+ (hadronic 
energy minus momentum), and Sh (hadronic invariant mass squared). The corresponding 



normalized distributions are given by 



mB-2Ei 



1 dT 



dCj 



r dEi 
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1 dV 
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oo 



r dsH 



1 dT 



TUB 



1 




(18) 



where gip) = g{A — u), and in the latter case 



F{uj, r) = 2r2(3 - 2r) S{lj) + 4r(l - r) £7ra j| g{td) . 



(19) 



Note that the contribution to the charged-lepton spectrum is much larger than that in the 
other two cases. We have studied the numerical impact of the four-quark contributions to 
these spectra using the model functions in ()14|) with A = 0.63 GeV and the default choice 
for the leading-order shape function from [0] . We take J'b = 200 MeV for the S-meson decay 
constant and as = 0.3 for the strong coupling at the intermediate scale. Even without assuming 
a significant color suppression the effects are very small. For e = 1, sizable distortions of the 
P+ and Sh spectra occur below P+ ~ 0.5 GeV and sh ~ 1-5 GeV^, while those of the charged- 
lepton energy spectrum are located above Ei ^ 2.2 GeV. For more realistic values \e\ <^ 1, 
the effects are even smaller. Varying the input parameters within reasonable limits does not 
change this conclusion. The four-quark contributions are proportional to the combination 
eTra^/^/A^ times a dimensionless function of the ratio cD/A. Assuming = (200 ± 30) MeV 
and A = (0.63 ±0.07) GeV, we have evraJl/A^ = (0.15 ± 0.07)6 GeV-\ where the parametric 
uncertainty can be absorbed into our ignorance about the color-suppression factor e. 

Figure |2l shows results for the partial rate fractions obtained by integrating the spectra 
in (fT^ over Ei > Eq = {uib — Ae)/2, P+ < Ap, and < sq = ttibAm, where the 
quantities Aj are chosen such that in all three cases the charm background starts at Aj = 
mj^/mB ~ 0.66 GeV. The effects of the four-quark shape functions are included under the 
extreme assumption of no color suppression {e = 1). In reality, we expect their effects to 
be significantly smaller. If we consider the partial rate fractions obtained by integrating the 
spectra over the phase space not contaminated by charm background, we find that the largest 
effects in the two models lead to 



Fe, = (6.5 + lAe)% , Fp^ = (60.9 - 2.3e)% , F,^ = (80.6 - 0.7e)% . (20) 



In the latter two cases these corrections are negligible even if the color-suppression factor e is 
not particularly small. 

For comparison, we note that the authors of ■SJ estimate (using naive dimensional analysis, 
and ignoring color suppression) that the impact of four-quark shape functions on the charged- 
lepton energy fraction could be as large as 180% of the leading term. Even for e = 1 this 
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Figure 2: Model predictions for the partial B — > X„/~z/ rate fractions with 
cuts on El, P+, and sh, respectively. The dashed line shows the leading-order 
result, while the two solid lines include the contributions from four-quark shape 
functions, evaluated setting e = 1. 

would be an order of magnitude bigger than our estimate. The main reason is that these 
authors include an "enhancement factor" g"^ = Airag ~ 4 in their estimate of the power- 
suppressed effects, as opposed to, say, nas ~ 1. Clearly, dimensional analysis cannot tell the 
difference between these two factors. The precise numerical coefficients multiplying the four- 
quark shape-function contributions to a given observable can only be determined by explicit 
calculation, as done in (|18|). Note that following the reasoning of fSJ one would associate an 
enhancement factor g"^ = Airag ~ 12 with soft-gluon exchange (for which ~ 1), which would 
be wrong, since the corresponding powers of gs are already included in what is conventionally 
called Aqcd- 

In summary, we have presented a simple, but well-motivated model for the four-quark 
subleading shape functions /^(u;) and fvij^), which contribute at order Aqcd/"^?; to inclusive 
5-decay spectra in the endpoint region. We have shown that in the vacuum-insertion approx- 
imation non-zero contributions only arise from four-quark operators for which the light quark 
flavor matches that of the 5-meson spectator quark. Vacuum pair production of the light 
quark pair, though not suppressed on general grounds, does not contribute due to the Dirac 
structure of the relevant operators. The result for the subleading shape functions /„ and is 
given by a color-suppression factor times a double convolution integral over a product of two 
5-meson light-cone distribution amplitudes. Using two simple models for this function, we 
have obtained explicit forms for the subleading shape functions, which are compatible with all 
known constraints from analyticity and moment relations. The corresponding impact on the 
decay distributions in semileptonic B — »• Xul~v decays is found to be negligible. While the 
results reported here are admittedly model dependent, we believe that they support our earlier 
claim [9J that these four-quark contributions are likely to be smaller than other subleading 
shape-function effects. 
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